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1 Introduction and summary 



There is a general understanding in the hterature that the LEP measurements of the gauge cou- 
phngs at the scale Mz and a remarkable convergence of the couplings at the scale around 10^^ GeV 
is a hint to a one step supersymmetric unified theory. There are however several physical arguments 
suggesting that in a larger unified theory like SO (10) there may be an intermediate scale cor- 
responding to a left-right gauge symmetry breaking somewhere around 10^^ to 10^^ GeV based 
on neutrino physics [Q, Q as well as strong CP problem [Q. However, such a scale must not affect 
the gauge coupling unification constraints. Recently a number of studies have been performed to 
examine this question ||, ^, |T^, |T^]. In this paper we present a systematic two-loop analysis of a 
class of intermediate-scale supersymmetric SO(IO) models varying the value of the strong coupling 
constant a^{mz) in the experimentally allowed range. To begin with we calculate the predictions of 
(i) unification scale Mx, (ii) intermediate scale of B-L symmetry breaking M/ (iii) the unification 
coupling adMx) and (iv) the proton lifetime in a two- loop renormalization group study. 

The large value of the top quark mass measured by the CDF and DO collaborations at FERMI- 



LAB [|T2| is suggestive of a large top Yukawa coupling. Keeping this in mind we include the effect 
of all the third generation Yukawa couplings in the running of the gauge couplings. We note that 
in the large tan j3 region the effect of the bottom quark and the tau lepton Yukawa couplings also 
affect the gauge coupling evolution non-trivially, and one-step unification prefers a lower value of 
a^{Mz) = 0.120 instead of 0.124 obtained in the low tan/5 region. This effect also exists in the 
intermediate scale models. Later in this paper we focus our attention to the masses of fermions. 
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We calculate the (v) pole mass of the top quark when the top quark Yukawa coupling is at the 
quasi-infrared fixed point at the scale Mx- The predictions of (vi) the mass ratio mh{mh) /mT-{mT-) 
also emerges from our analysis. The dependence of these predictions are displayed graphically 
and they are mostly consistent with the observed values. 

The ratio of VEVs tan f3 = (H^) / {H'^) is an unknown parameter of the supersymmetric models. 
However it has been noted that for very small values of tan/3 the top quark Yukawa coupling 
becomes very large before the scale Mx and consequently it is possible to give a lower bound on 
tan/3 requiring the perturbative unitarity of the top Yukawa coupling up to the scale Mx- In this 
paper we make an analysis on the (vii) lower bound on tan/3 in such intermediate scale models. 
Taking the top quark Yukawa coupling to be arbitrarily large at the scale Mx we derive an upper 
bound on it at the top scale [ht{'mty\ and then we convert this upper bound to a lower bound on 
the unknown parameter tan j3 for given values of the top quark mass. 

We know that is difficult to understand the origin of a small neutrino-mass in the minimal 
supersymmetric standard model (MSSM) or in the minimal SU(5) GUT. In the presence of a right 
handed singlet neutrino the VEV of the SU{2)l doublet Higgs scalar will provide a Dirac mass to 
the neutrino which comparable to the masses of the other fermions. On the other hand we also 
know that in an intermediate scale model such as SO(IO) GUT, it is possible to relate the mass of 
the left-handed neutrino to the inverse of the large lepton number violating Majorana mass of the 
right-handed neutrino {Mn ~ Mj) by a see-saw mechanism |Q. We study the (viii) mass of the 
left-handed neutrino (of the third generation) in these models and calculate predictions for the tau 
neutrino mass by inputing of the values of the intermediate scales and various Yukawa couplings 
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calculated from the first part of the paper. All our results can be compared with the one-step 
unification scenario in the limit M/ = Mx- 



We will consider the value of a^[Mz) |13l in the range 0.110 to 0.130, as the LEP measurements 



give a2,{Mz) = 0.124±0.006 [|T4|, which is typically 2 to 3 standard deviations larger than the value 



coming from J/\E'-decay which is az{Mz) = 0.108 ±0.010 [15|, from lattice calculations which gives 



the value a^{Mz) = 0.110 ±0.006 |16], or from deep inelastic scattering experiments |17|which give 



a^{Mz) = 0.112 ± 0.002 ± 0.004. 

This paper is organized as follows. In section (II) we motivate our models by a brief summary 
of an one-loop analysis. In section (III) we give the two-loop analysis of the mass scales. In section 
(IV) we give the predictions of top and the bottom quark masses and derive lower bound on tan f3. 
In section (V) we discuss the predictions of tau neutrino mass from see-saw mechanism and in 
section (VI) we conclude. The relevant formulae have been summarized in the appendix. 



2 Models and scalar contents: A case study 



Before doing the rigorous analysis it is helpful to do a one-loop case-study and fix the scalar contents 
of various models to be considered afterwards. It is important to note for our purpose that an unified 
model with single intermediate scale is as predictive as a model with grand desert in the following 
sense. Three low energy observables in the gauge sector namely ai(M^), a2{Mz) and a^^Mz) can 
determine the three unknown parameters namely Mx, Mj and adMx)- Given ai{Mz) and 0:2 (M^) 
sufficiently precisely, there exists a unique value of 03 (M^) in the experimentally allowed region for 
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which one gets Mx = Mi. Indeed in a model with more general scalar and fermion content one 
gets Mx 7^ Ml, and thus, here we seek the models in which Mi < Mx- 

In one loop approximation this unique value, for which Mx = Mi, is a^{Mz) ~ 0.1144. It is 
noted by Lee and Mohapatra that it is possible to get an intermediate scale M/ < Mx in a 
number of models if 7^ 0.1144. In this paper one of our objectives is to find this unique value of 
03 in two-loop approximation and derive the ranges for the intermediate scale for the allowed range 
of a^{Mz). 

Now let us write down the 1-loop renormalization group equation (RGE) of the three couplings 
introducing a general intermediate scale M/ between Mz and Mx- We have used hi to denote the 
beta function coefficients below the intermediate scale and h[ to denote them above the intermediate 
scale. The relations are. 



A combination 5 can be written [TH] in which hi get eliminated but h[ survive, as. 



5 = 7a^\Mz) - 12a^\Mz) + ha^\Mz)- (2) 
Eqn.(0) and Eqn.(^) together leads to, 

5 = -{7b', - 126', + 56;) In ^ = - In ^. (3) 
27r^ ^ ^ Ml 27T Ml ^ ^ 

Now, if ai(Mz) = 0.01696, 0^2 (M^) = 0.3371 to get one-step unification in one loop we need 
a^lMz) = 0.1144 which gives 6 = from Eqn.(0). However, experimentally a^lMz) can be in the 
range 0.11 to 0.13. Hence we obtain an allowed range in A which is crucial to get an intermediate 



scale, 

- 2-K Mi - ^ ' 

It may be noticed that when A is small we can achieve a large value for In leading to an 
intermediate scale well-below the unification scale. We will return to the inequality ^ in a moment. 
Now let us consider the symmetry breaking pattern, 

50(10) 5f/(3)e X 5t/(2)i X 5f/(2)fi X f/(l)(B_i), 

^ SU{?>),xSU{2)lxU{1)y, 

^ 5f/(3)e X f/(l)e„. 

We restrict the type of Higgs representations above the intermediate scale Mi by requiring that 



the supersymmetric SO(IO) GUT emerges from an underlying superstring theory |19|. Restricting 
ourselves to only those Higgs scalars which can arise from simple superstring models with Kac- 
Moody levels one or two, we can have a restricted number of solutions of Eqn.(^. The various 
models can be characterized by a set of five integers {nL,nii,nH,nc,nd), where nc refers to the 
number of (8,1,1,0), uh means the number of (1,2,2,0) fields and ul and ur means the number 
of (1,2,1,1) + (1,2,1,-1) and (1, 1, 2, 1) + (1, 1,2,-1) fields under the intermediate symmetry gauge group 
and rid refers to the number of pairs of Higgs doublets below the scale M/. To further restrict the 
number of models we have the following assumptions. 

• The theory below the intermediate scale in strictly the Minimal Supersymmetric Standard Model 

(MSSM). We will also consider an example where the theory below M/ has four Higgs doublet^ 
In this case the combination in the Eqn. (||) is no longer orthogonal to the beta function coefficients be- 



instead of two in the standard model. 



• It should be possible to get an intermediate scale which is at least one to two orders of magnitude 
below the unification scale ( In = 2.3 — 4.6). 



Using the beta function coefficients listed in the appendix Eqn. (||) gives, 

- 10 < (-9 + 21ne - 9nH + Qur - dni) < 6. (5) 

The models satisfying the above inequality are tabulated in Table |I| which are numbered in the 
decreasing order of minimality. We have included one scenario namely model VIII which has a 
color octet above the scale M/. This scenario has A = and thus it can accommodate a large 
splitting between Mj and Mx- In scenario IX the theory below the scale Mj has four Higgs 
doublets {ud = 2). 



Model 


riL 


riR 




ric 


rid 


A 


I 





2 


1 





1 


-6 


II 





3 


1 





1 





III 





4 


1 





1 


6 


IV 





3 


2 





1 


-9 


V 





4 


2 





1 


-3 


VI 





5 


2 





1 


3 


VII 


1 


5 


2 





1 


-6 


VIII 


1 


1 


1 


1 


1 





IX 





3 


2 


1 


2 





Table 1: The minimal models which satisfy the condition in Eqn. (5). When the quantity A is 
positive (negative) the model gives rise to an intermediate scale for the lower (higher) values of 
a^lMz) than the one step unification case for which a^lMz) = 0.1144 at the one-loop level. In the 
case A = the intermediate scale is unconstrained at the one-loop level. 



low Mi. However it is straight-forward to derive an equivalent combination for this case. Interesting low-energy 
phenomenology of such a model can be found in Ref . [pO| . 
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3 Two loop analysis of mass scales, unified coupling and 
proton decay 



In the two-loop case the gauge couphngs evolve according to the equation, 

where the Yukawa couplings can be defined by the superpotential invariant under the intermediate 
symmetry, as, 

= hQ,Q'^T2(l)lQ'' + hQ,Q^T2<l)2Q'' + hL,L'^T2(f)lL'' + hL,L'^T2(f)2L'' + h.C., (7) 

where we have denoted the quarks and leptons by the obvious notation Q, Q'^ and L, L'^. At the 
GUT scale, we have hQ,{Mx) = hL,{Mx) - hi{Mx) and hq^^Mx) = Hl^^Mx) = h2{Mx). (pi,2 are 
two Higgs bidoublets embedded in 10i,2 scalars of SO(IO). We have used the notation, 

t = \n^i , a = - , Y = -, (8) 

where g and h refers to the gauge and Yukawa couplings respectively. The coefficients are given 
in detail in the appendix. When Yi{Mx) and Y2{Mx) take large values of order 0(1) the Yukawa 
couplings affect the running of the gauge couplings most. In this case Yt and Y;, have almost equal 
values at the low energy scale. Such a situation corresponds to large tan /3. On the other hand if 
Y2{Mx) « Yi{Mx) at the scale Mx the top-bottom mass hierarchy is due to Yukawa coupling 
hierarchy and consequently tan /3 is order 0(1). We will refer to this case as the small tan P scenario. 
In models I, II, III and VIII tan /? is always large. 

First, we present the results of the unification scale. The variation of the unification scales with 
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that of as{Mz) have been plotted in Figure (la). The predictions of the intermediate scale have 
been plotted in Figure (lb) and that of the unification coupling adMx) have been plotted in Figure 
(Ic) for various models. The grand desert case can be recovered from the meeting point of all the 
curves in Figure (lb) that is when the intermediate scale is equal to the GUT scale. In the low 
tan (3 region the meeting occurs at the value as{Mz) = 0.124 and in the high tan (3 region it occurs 
at the valueQ a^lMz) = 0.120. For all other values of a^lMz) it is possible to have an intermediate 
scale. In model VI the unification scale becomes low in the low as region; for model V the same 
thing happens for high as{Mz) region. As the dimension five proton decay can be suppressed in 



the SO(IO) models by some additional mechanism |22|, we plot the dominant dimension six decay 
mode in Figure (Id). The proton decay rate has been calculated using the formula below in which 
we have taken the mass of the heavy gauge bosons as Mx- 

r(p ^ ^^7) = i^^^? 11 + (1 + ly^^ri n (9) 

The values of parameters used in Eqn. (^, have been listed in the appendix. The variation of 
decay mode vis-a-vis the variation of a^lMz) comes not only from the variation of Mx and ac 



but also due to a (mild) renormalization of the proton decay operator ||2^, ^ up to the scale of 
1 GeV which has been parameterized in the factor A^. When a^lMz) is above 0.130, as{lGeV) 
becomes very large and the perturbative calculation is not dependable. In our estimation of Aji we 
have made taken into account the renormalization of the operator due to SU{3) color and SU{2)i 
interactions only. This is a reasonable approximation as in most cases M/ is large enough and the 

SU{2)r renormalization effects exist only beyond M/. In the one-step unification case our value of 

For moderate values of tan/3 unification occurs for values of aj,{Mz) within 0.120 and 0.124. Low energy 
threshold corrections pl| assuming universality of gaugino masses will shift this whole range upwards. 



the parameter Ar agrees with that of Ref. 



The hfetime of proton has been plotted in Figure (4). Model V predicts large proton decay 
rates (in the borderline of experimental limit) for the values of as{Mz) beyond 0.123. Similarly 
model VI predicts equally large rate of proton decay for lower values of a^lMz) below 0.114. The 
models except I, VIII and IX predicts a lower proton lifetime than the minimal one step unification 
scenario when there is an intermediate scale. This is encouraging from the forthcoming proton decay 
experiments where these models can be probed. In fact if the experiments at SuperKamiokande 
excludes the minimal unification scenario it will also exclude all the present models except model 
I, VIII and IX. 

Before concluding this section let us make a small remark about the curvatures visible in the 
plots near the meeting point of the models. These curvatures are more pronounced in models V and 
VII. Let us consider Figure (lb) for example. Near the GUT scale the Yukawa couplings are large 
and they fall quickly below the GUT scale. The Yukawa couplings tend to pull the individual lines 
towards the low region wheres in the high models (like V and VII) the gauge interactions 
have exactly the reverse effect. This causes the curvature in the graphs near the GUT scale which 
is purely a two-loop effect. 

4 Fermion masses and Yukawa couplings 



The simultaneous evolution of gauge and Yukawa couplings enables us to calculate the following 
Dirac type masses of quarks and leptons, in respective scales denoted in the parenthesis, by the 
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following relations, 



mt{mt) = ^^4:nYt{mt)-j=sinp, (10) 



mbinif,) = ATiYbimt) -j= cos P rjb, (11) 
I V 

mr{mr) = j4:7rYr{mt)—^ COS P 1]^-. (12) 

V 2 

m,^{Mj) = ^A'KYLi{Mi)^sm(3. (13) 
The factors r/b and rj^- takes into account the running of the masses of mb and uij. to their respective 



scales starting from m^. We have used the formula of reference to calculate ?7b and r]^ taking into 
account the three loop QCD effects and one loop QED effects. We have got rj^ = 1.017 assuming 
a~^{Mz) = 127.9. The values of rjb are plotted in the appendix. We start from a pair of values 
Yi(Mx) and Y2{Mx)- Using the RGE for the Yukawa couplings given in the appendix we calculate 
the values of the couplings Yt{mt), Yb{mt) and Yr{mt). At the scale Mj the boundary conditions 
for the Yukawa couplings are, 

Fqi (Mj) = Yt (Mj) , Yq2 (Mj) = Yb (Mj) and Yl2 (Mj) = Y, (Mj) when hh = 2, (14) 
Yq (Mj) = Yt (Mi) = Yb (Mj) , and Yl (Mi) = Y, (M,) when uh = 1. (15) 



Using Equation Eqn. (^2|), which has the least experimental error we calculate the value of cos/3. 
Once the value of cos (3 is known the predictions for mt{mt) and mb{mb) follows from Eqn. (|1^) and 
Eqn. (|ri|). The running mass of the top quark has been calculated by iterative procedure using the 
condition 

mtimt) = mt. (16) 
11 



The pole mass |]23| has been calculated from the running mass in Eqn. ( [I6| ) using the equation 

pole . Nr., 4^3 (mj 

m^^ = mt(mj) [1 H J, (17) 

OTT 

for each value of a^^Mz). The predictions of mf°'*^ is plotted in Figure (2a). The value of Rh{mt) = 
are plotted in Figure (2c), and when the prediction of ■mb(mt) is extrapolated to the mass 
scale of the bottom quark we get Figure (2d) using mT-(m^) = 1.777 GeV. A few comments are in 
order. 

• When we decrease sin/5(tan/3), the reduction in the bottom quark Yukawa effect increases the 
value of htlrrit). This is visible in the plot (2a). However, once the value of Y2{Mx) is smaller than 
10""^, the effect of the bottom quark Yukawa coupling on Yt virtually ceases to exist. At this stage 
a further reduction in sin/3(tan/3) causes a large reduction in the predicted value of the top mass. 
And in fact for Y2{Mx) = 10~^ the predicted top mass falls down to 161 GeV. 

• The experimentally allowed range of mf,(mb) is between 4.10 and 4.50 GeV flSl. However masses 



as high as 5.2 GeV are also considered in the literature. The predictions put strong constraints on 
the model. Especially for the low aslMz) case, high tan/? regime is disfavored. In the low tan/? 
scenario the predictions shift towards higher a^iMz) which leads to a bigger rjh [see Fig. (7)] and 
hence a larger predicted value of mbirnb). 

• The lower bounds on tan/3(sin/3) follows directly from an upper bound onQ Yt{mt). We calculate 

the bound by the following procedure. Taking the value of Yt{Mx) to be arbitrarily large we 

calculate the upper bounds for Yt{mt) for various Y2{Mx)- The absolute upper bound is converted 
^It is well-known that the Yukawa evolution equations have an infrared quasi-fixed point structure [ p7[ . 
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into a lower bound for tan/? by Eqn. (p^ and Eqn. (p!7|). These bounds have been plotted in 
Fig. (2b). 



5 See-saw mechanism and neutrino mass 



If the neutrino is a Majorana particle it can have a lepton number violating Majorana mass term. 
SO(IO) GUT has a natural mechanism to generate a large Majorana mass of the right handed 
neutrino. In this case the neutrino mass matrix has the form, 



/ m'^\ 



M 



\m'^ M J 

where m'^ is a 3 x 3 Dirac mass matrix and M is a 3 x 3 Majorana mass matrix in general. When the 
matrix Ai is diagonalized we get two eigenvalues of orders M and [mf'Y /M . The latter eigenvalue 
can in principle explain the smallness of the mass of the left-handed neutrino when M is large. 

In our models the intermediate B — L symmetry is broken by the Higgs scalars 16 + 16 fields 
of SO(IO). We will consider two different scenarios by which Majorana mass of the right handed 
neutrino can be generated. 

(a) Using a higher dimensional operator of the form -^IQplQplQulQH written in terms of SO(IO) 
representations. The subscripts F and H mean fermions and scalars respectively. When IQh gets a 
VEV a large Majorana mass of the order h v\/Mx is generated. 

(b) Introduction of additional singlets to have a generalized see-saw mechanism In this case 
the mass is effectively generated from a 3 x 3 see-saw matrix. This case is interesting as an effective 
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2x2 see-saw matrix can be recovered purely from renormalizable interactions when the singlet 
scalars get their VEVs. 



The prediction of the mass of the left-handed neutrino in the case (a) is well-known. We recast 
it in terms of the parameters being evaluated by RGEs in our cases (we have used M/ = vr) , as, 



h m^{Mj) = z:t1T7 — = 87r^a2/?(M/)- 



X 



Mf/M_ 



(19) 



X 



In the case (b) |28[, the 3x3 mass matrix in the basis ((//„, Na, Sa), where a=l,2,3 and Sa are 



the singlets) is, 



M 



( /lit;" \ 



hiv"" h'vR 



(20) 



V h'vR Ms ) 

where is the mass of the singletR. In this case the small mass of the left-handed neutrino arises 



purely from the renormalizable interactions, as. 



h' m^{Mi) 



Yli{Mj) [t>nan=^/j/(l + tan^/?)] 
' MfjMs ' 



(21) 



In Eqn. (^TD we will consider Mg = Mj. When Mg = Mx we recover Eqn. (|T^) from Eqn. (^T]). In 



both the cases we will assume, 



247V2 GeV'^ 



(22) 



In Figure (3a) and Figure (3b) we have plotted the left handed neutrino masses in scenarios 

(a) and (b) modulo the unknown Yukawa couplings h and h' in various models as a function of 

as{Mz)- These numbers will get renormalized when they are extrapolated up to the scale of 1 GeV 
^For the detailed superpotential see RefS 
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29| , |30| , pl[] . In Ref |31| it has been noted that the Yukawa couphngs affect the extrapolation of the 
see-saw formula and the tau neutrino mass increases. This extrapolation will change the mass of 
the left handed neutrino by a factor of the order of unity. 

We know that a tau neutrino mass of the order of a few electron volts is preferable if neutrino 
is to be a candidate for the Hot Dark Matter (HDM). We see that in scenario (a) a tau neutrino 
in the range of 1-10 eV can be achieved in all the models depending on the value of 0:3 (M^). On 
the other hand, scenario (b) can predict a tau neutrino mass in the 1-10 eV range for models V, II 
and VI. However, model VI and V have potential problems with proton decay around the values of 
a^i^Mz) needed to produce a correct order of the mass of the u^-, whereas model II predicts a lower 
value for the bottom quark mass. 

6 Conclusions 



We have considered a class of minimal models for which an intermediate scale in a SO (10) GUT 
can be achieved. We have found that a small number A can be defined, the sign of which can tell 
whether the unification occurs for a low value of 0:3 or a high value of compared to the one-step 
unification. After listing the minimal models in which an intermediate scale can be generated, we 
have done a two-loop analysis of the gauge couplings to find the unification scale, the intermediate 
scale and the unification gauge coupling. These predictions are in general functions of the input 
a3{Mz). We have plotted the predictions in Figure 1. The predictions for one-step unification can 
be recovered from the meeting point of all the branches in these figures for which Mj = Mx- The 
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gauge boson mediated proton decay rates have been calculated after renormalizing the proton decay 
operator to the scale of IGeV. These predictions are also plotted in the Figure 1. In some models 
the predicted proton life-time is in the borderline of experimental exclusion limits. These models 
can be tested in the SuperKamiokande experiments soon. 

We have also done a parallel two-loop analysis of the Yukawa sector of these models. The 
combined running of the gauge and Yukawa couplings enable us to calculate the predictions of the 
top and bottom quark masses using the mass of the tau lepton m^(mT-) as an input. The predictions 
of the top mass is always in the range given by the CDF and DO collaborations. The bottom quark 



mass is not always in the range 4.10-4.50 GeV as quoted in the review of particle properties [|T3 . 
Thus predictions of the bottom quark mass is a good indicator by which one could compare these 
models relative to each other. Running of the Yukawa coupling also gives a lower bound on the 
parameter tan/3 from the requirement that the top quark Yukawa coupling remains perturbative 
up to the scale of unification. 

The determination of the intermediate scale and the running of the Yukawa couplings enable us 
to calculate the predictions of the neutrino masses via see-saw mechanism. The Dirac type entry in 
the see-saw matrix can be fixed from the running of the neutrino Yukawa coupling. The predictions 
of the tau neutrino for various values of 0:3 have been plotted. A mass of the tau neutrino in the 
range 1-10 electron volts can be achieved in various models. It can be noted from the plot of the 
neutrino masses that in the one-step unification model the mass of the tau neutrino is in the range 
of 10~^ electron volts only, making it less attractive as a dark matter candidate. 
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The scenario VIII can be a good candidate for string unification with the unification scale near 
the plank scale and the unification coupling larger than the one-step unification case. 



Acknowledgments 

We thank R. N. Mohapatra and G. Senjanovic for useful comments. We also thank J. Kubo for 
sending us a numerical subroutine which we have used in our analysis. 

7 Appendix 

7.1 Beta function coefficients for the Yukawa couplings 



It is easy to write down the evolution equations for the Yukawa couplings defined by ^JiTiYijk (j)i4>j4>k, 
as, 

"^^''^ -2Yi,k bi + lj + lk]. (23) 



dt 

The scale t is defined as t = In fi. In the two loop approximation we can write, 

7i = 7f^ + 7f- (24) 

The superscripts in parenthesis denote the order of perturbation theory. The anomalous dimensions 
for the superfields below the scale Mj are given below. In the one-loop anomalous dimensions are, 

- (26) 
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% - 


1 

47r 


[lYt - -as - J^«iJ' 


1q - 


1 

47r 


8 3 
[Ft + Ffe - -as - -a2 


ts; = 


1 




t£ = 


1 

47r 


3 3 , 
Wt- -a2- —ai]. 



1 

30^ 



The two-loop anomalous dimensions are; 

1 



(2) 
iL = 



'E 



7g - T^[-^.(7?+7f)-3n(7g^+7f) 



(28) 
(29) 
(30) 
(31) 



T^[->^^(7ff+7S) 

-(^«2 + |^«i)7? + ^&2«2 + (32) 
^[-2^(71^7^) 

-^«i7g^ + ^M?], (33) 

^[-2n(7S^ + 7g;) 

8 2 8 2 

,8 8 , 8 2 8 21 

-(3^3- Y^ai)7^ +3M3 + Y^&iai], (35) 

= J^[-l'.(#+7S)-Vi(7g'+7g) 

-(^a3 + ^Q;2 + ^fciQ;i)7Q^ + ^^'s^s + ^^2^2 + ^^i^^il' 



(36) 
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;^a2 + ^ai)7Q^ + ^&2a2 + ^&iai], (37) 



;^Q;2 + ^ai)7S + ^^2a2 + ^6iQ;2]. (38) 



Above the scale Mj the one-loop anomalous dimensions are, 



= — [2yQi + 2rQ2--ac- 2«L- (39) 

7gi = i^[2FQi + 2FQ2-^ac-^aK-^«B-L], (40) 

7?) = i_[2r^^ + 2F^2-^«L-^«s-L], (41) 

7^ = ^[2lLi + 2F^2-^ai?-^as-L], (42) 

7if = ^[3% + r^,-^a^-^a«]. (43) 



The two-loop anomalous dimensions are, 



- Y^[-2V(7i;^+7S)-2yg2(7£^+7gi; 



3 1 ^ , 8, ^2 3, ^2 1 T ^2 



[-ac + -ftL + — aB_L)7y + -ha^ + -6^^^ + — ^B-Ltts-J, (44) 



16 TT 



,8 3 1 \ (1) ^ T 1 7 2 2n /^^\ 

{-ac + -an + —aB-Lhqi + -ha^ + -bROiR + —hB-LOiB-Ll (45) 



7r = Y^[-2F,i(7;r+7icO-2F^2(7;,^+7iV) 

-(^"L + ^(^b-l)1l^ + ^^L^i + ^bB-Lal_L\, (46) 
tS' = J^l-2r..(7;?+7i")-2y.,(7£>+7i") 
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-(^"fl + ^«b-l)7lc^ + ^bnal + ^6B_La|_i], (47) 



7i? = Y^h^.a7in7S)-3r,.(7? + 7g) 



|«L + ^«fi)7i? + ^hal + (48) 



7.2 Beta function coefficients for the gauge couplings 

The two loop RGE for the gauge couphngs are given in general by Eqn (^. The coefficients 6j, bij 
and depend on the scale range considered, and are listed below: 



(A) For Mz < < Mj the coefficients are: 

/ 2nf + l 



—6 + 2nf + Ud 
\ -9 + 2n/ 



1,2,3, 



(49) 



irid 



5^/ + s^-i 



V 5 



-24 + 14n/ + 7nrf 
3n/ 



15^/ 

-54 + f 



1,2,3, 



(50) 



26 


14 


18 


5 


5 


5 


6 


6 


2 


4 


4 






t,b,T. 



(51) 



(B) For Mi <fi< Mu- 
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/ 2nj + |K + n^) \ 
-6 + 2nf + nH + nL 
-6 + 2nf + nH + nji 
\ -9 + 2nf + 3nc J 



B — 3 . 



(52) 



\ 



3n/ + |ni 
-24 + Uuf + 7{nH + til) 
3nH 
3nf 



3nf + |nR 
3nH 

-24 + 14n/ + 7{nH + nR) 
3nf 



8nf 

-54 + ^n/ + 34nc / 



i, j = B — L, L, R, 3 . 



(53) 



/ 


2 


2 


6 


M 




12 


12 


4 


4 




12 


12 


4 


4 


I 


8 


8 





/ 



B-L, L, R,3, k = Ql, Q2, LI, L2 . 



(54) 



7.3 Proton decay parameters 



The proton decay formula is given by Eqn 



the parameters 33 : 



: where we have assumed the following values for 



D = 0.81 ; F = 0.44 ; = 139 MeV 

rriN = 0.938 GeV ; |Kd| = 0.975 ; a = 0.03GeV^. 

The value of a is taken in such a way that proton decay is maximal. In fact, its value has a range 
from 0.03 to 0.003 GeV^. The value of the renormalization factor varies from 2.7 to 4.1 The 
values are plotted in as can be seen in Fig. (4). 
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7.4 The range of rji^ 



Here we give the values of rjb used to extrapolate mb{mt) to mb{mb), we have calculated its 



variation with 03 (M^) for the different models studied in the text in Fig. (5). We have taken 



ct2{Mz) = 0.3371 and aiy(M^) = 0.01696. Our calculation agrees with that of Ref. ||26| quite well. 
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Figure 1: Predictions for (a) Unification scale, (b) Intermediate scale, (c) Unification gauge cou- 
pling, and (d) Proton life-time, for the models listed in Table I. Solid lines denote high tan/3 
{Yi{Mx) = Y2{Mx) = 1; dashed lines denote the low tan/5 regime {Yi{Mx) = 1, Y2{Mx) = 10"^). 
In Figure (d) the dotted line is the experimental limit Tp = 5.5 x 10^^ yr. 
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Figure 2: Predictions of (a) pole mass of the top quark, (b) lower bound on tan/3, (c) R^t at m^, 
(d) running mass of the b quark. Solid lines and dashed lines are as in Fig. (1). 
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Figure 3: Predictions of the left-handed neutrino mass by see-saw mechanism by the two scenarios 
(a) and (b). Sohd and dashed hues are as before. 




Figure 4: Variation of with respect to a^lMz) in different models. Solid and dashed lines are 
as before. 
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Figure 5: Variation of 77^ with respect to a-^iMz) in different models. Solid and dashed lines are 
as before. 
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